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OJ ■ Abstract 

^ ' We discuss various methods and their effectiveness for solving linear 

equations over finitely generated abelian groups. More precisely, if ip: G — > 
I/"") ' H is a homomorphism of finitely generated abelian groups and b 6 H , we 

discuss various algorithms for checking whether or not b 6 imcp holds and 
if so, for computing a pre-image of b in G together with the kernel of ip. 

o 

Solving linear equations over finitely generated abelian groups is an important 
tool for computing H 2 (G,A) of a given G- module A. Therefore algorithms 
for solving linear equations are important in extension theory of groups, see 
Section 8.7 of [6]. For extending the Small Groups Library, see e.g. [1], we need 
t> ' effective algorithms especially for this problem. In this note, we discuss various 

methods and we propose an effective algorithm which performs well on a broad 

■»h- . range of randomly chosen linear equations. 

(N 

q ■ 2 Preliminaries 

O" 

Let G and H be finitely generated abelian groups and let ip: G — > H be a 
homomorphism. Further let b G H be given. A linear equation over G asks 
whether or not b £ map holds and, if so, for computing a pre-image of b in G 
^ ■ together with the kernel of ip. 

3 Standard methods for polycyclic groups 

Since the groups G and H are finitely generated and abelian, both are polycyclic. 
Therefore the methods for homomorphisms of polycyclic groups described in [2J 
Section 3.5.2] apply. We briefly summarize these methods here for completeness. 

The homomorphism (p:G^H naturally corresponds to the subgroup Uu, = 
{{d^ i9) I 9 £ G} of the direct product H xG. Given polycyclic sequences Q and 
T-i of G and H respectively, the algorithms in [2] allow to compute an induced 
polycyclic sequence for U^ with respect to 7~L x Q. This sequence allows to read 
off the kernel of <p and can be used to decide whether or not b E vunp holds. 



Lemma 1 (Eick, 2001) Let ip be a group homomorphism and hi be an in- 
duced poly cyclic sequence for U^ with respect to % x Q. Then U has the form 
((ui,ui), . . . , (u s ,u s )) with Ui € G and Ui € H. Let t be maximal with Ut ^ 1. 

(A) Then (u\, . . . ,Ut) is an induced poly cyclic sequence for the image of (f with 
respect to TL. 

(B) Then (itt+i, • • • ,Us) is an induced polycyclic sequence for the kernel of if 
with respect to Q. 

Proof. For a proof, we refer to [2j p. 36]. □ 

This method is implemented in the POLYCYCLIC package of GAP, see [3]. How- 
ever, for finitely generated abelian groups, we can do much better as we will 
show in the following. 

4 Solution by solving linear Diophantine equations 

In this section we describe a more effective algorithm relying on linear algebra 
only. Let G be a finitely generated abelian group. Then G decomposes into its 
torsion subgroup T(G) and a free abelian subgroup so that G = Z x T(G) holds. 
Further, the torsion subgroup T(G) decomposes into its p-Sylow subgroups so 
that we may identify 

g = i} e s y i pi (T(G)) e • • • e s y i Pfe (t(g)), 

where SyL (T(G)) denotes the pj-Sylow subgroup of T{G). Similarly, we can 
identity 

H = J}' Sy\ pi (T(H)) e • • • Sy\ pk (T(H)). 

Every homomorphism ip:G—>H induces a homomorphism ipi = <p\s y \ (T(G)) °f 
the j?j-Sylow subgroup. Assume that Syl pi (T(H)) = Z «« © • • • © Z a in . holds. 

Then, with respect to an independent generating set of G, the equation x v = b 
translates to equations of the form 

A q a s )^ = (3 (modfpr 1 ,...,^— ,0,...,0]) (1) 

where anx\ + • • • + ai n x n = bi (mod 0) denotes the Diophantine equation, 
the sub-matrix At is a block diagonal matrix describing the homomorphisms 
tpi , . . . , (fk , while A s incorporates the action of ip on the free abelian subgroup 

Z e ^G/T{G). 

For solving ([T]), we introduce free variables, one for each generator of the 
torsion subgroup T(G). Thereby we obtain a linear system of Diophantine 
equations of the form 

A t B D 



A s 

where D = diag(p" n , . . . ,p amn ™) is a diagonal matrix describing the modular 
equations of the form anx\ + • • • + aaX£ = bi (mod p ai i). 



4.1 Solving linear Diophantine equations 

The linear system of Diophantine equations in ([2]) can be solved using the Smith 
normal form; see [7]. Recall that, for every A E % m ^ n there are unimodular 
matrices L € SL(m,Z) and R € SL(n, Z) such that LAi? is a diagonal matrix 
D = diag(di, . . . , d r , 0, . . . , 0) whose non-zero entries di, . . . , d r satisfy di | <ij+i 
for every i < r = rk(yl). 

Therefore the linear system Ax = b of linear Diophantine equations in ([2]) 
is equivalent to the equations Dy = Lb and x = Ry. Write c = Lb. Then the 
system Dy = c has an integral solution if and only if di | c% holds for every 
1 < i < r and Cj = otherwise. If the equation Dy = c has an integral solution, 
then these are given by y = (c\/d\, . . . , c r /d r ,ti, . . . , t n _ r ) with t\, . . . , t n _ r G Z. 
Every solution x to the linear system of Diophantine equations Ax = b is then 
easily obtained asx = Ry. 

It is well known, 0, that computing the Smith normal form of an integral 
matrix is computationally hard due to the unavoidable growth of the interme- 
diate matrix entries. Even though there are improvements in special cases, see 
or Section 9.3] for an overview, we cannot expect the algorithm in this section 
to be practical in general. Especially for groups with a large torsion-subgroup, 
the problem grows significantly by forming the system of linear Diophantine 
equations in ([2]) . Therefore reducing to a minimal generating set for the torsion 
subgroup first might improve this algorithm. Note that especially finite groups 
are a problem here. 

We can use a different and more effective approach for finitely generated 
abelian groups which first considers the linear system of Diophantine equations 
arising from the torsion-free part in G and H. Afterwards, our approach solves 
the equations over the torsion subgroups with varying the right-hand-sides with 
respect to the solutions of the torsion- free part. 

5 Solving linear equations over finite abelian groups 

In the remainder we consider linear equations over finite abelian groups only. Let 
G and H be finite abelian groups. Suppose that pi,...,p n are prime numbers 
such that |G| = p^ 1 • • ■pff and \H\ = p^ 1 ■ ■ -p„ with ej, fi 6 Nu{0}. Since every 
homomorphism cp-.G^-H induces a homomorphism i^j: Syl p . (G) — >• Syl p .(H) of 
the pj-Sylow subgroups, we may restrict to the case that G and H are abelian 
p-groups. More precisely, if we identify 

G = Syl pi (G)x...xSyl pn (G) and H = Syl^ (G) x • • • x Syl^ (G) , 

then we can decompose ip into homomorphisms (/jj:SyL. (G) — > Syl p .(H) with 
vlsyi (G) = f%- This yields that ip decomposes as (fix- ■ -X(p n , where ipiX- ■ -X(p n 
acts diagonally by 

( gi ,...,g n r*-^=(gf\...,g^). 



Computing pre-images and the kernel of tp can be done independently for the 
Pi-Sylow subgroups. In particular, solving linear equations over finite abelian 
groups split into independent computations for the Sylows subgroups and hence, 
can easily be parallelized. 

6 Solving linear equations over abelian p-groups 

Let G and H be abelian p-groups and let tp: G —■ H be a homomorphism. We 
identify G = 7L V ^\ x • • • x Z p e„ and H = Z / x x • • • x Z„/ m with e± < . . . < e n 
and /i < . . . < f m . Further let {g\, . . . ,g n } and {h\, . . . , h m } be independent 
generating sets of G and H with \gi\ = p e ' and |/ij| = p-", respectively. Then 
we can represent tp by an m-by-n matrix A = (a>ij)i,j where (au, . . . ,a m j) = 
{pt\, . . . , a m ) whenever g^ = /i" 1 • • • /i^ m holds. 

Computing a pre-image of b = h^ ■ ■ ■ h^ £ H together with the kernel of tp 
is equivalent to solving the linear system of modular equations 

anxi + ••• + ainXn = h (modp-^ 1 ) 

:•-.:! ! (3) 

a m ixi + ••• + a mn x n = b m (mod p fm ) 

We discuss the following methods for solving (|3|): 

(A) solving a linear system of Diophantine equations, 

(B) using the method for polycyclic groups described in [2], 

(C) using a modular analog of Smith normal form in the case that G and H 
are both homocyclic, 

(D) using the method in (C) and lifting solutions recursively, 

(E) lifting a solution over ¥ p recursively with the Hensel lemma. 

In Section ?? we show the application of these algorithms to various randomly 
chosen equations. 

6.1 Solving linear equations over homocyclic p-groups 

Let G and H be homocyclic p-groups; that is, G = Z i x • • • x Z « for some k 
copies of Z e . We generalize the algorithm of Section 0] to solve linear equations 
over homocyclic groups. 

Let tp: G — >• H be a homomorphism and let g G H be given. Then the 
endomorphism tp is represented by an m-by-n matrix A = (aij)ij while g G H 
translates to its corresponding exponent vector (bi, . . . ,b m ). For computing all 
solutions to the system Ax = b (mod p e ), we use the following modular analog 
of the Smith normal form: 

Lemma 2 Let A € (Z e) mxn be given. Then there exist matrices L G (Z i) mxm 
and R G (Z p t) nxn , which are invertible modulo p^, so that LAR is a diagonal 
matrix D = diag(di, . . . , djt, 0, . . . , 0) whose non-zero entries d±, . . . , d/. satisfy 
di | dj+i for each 1 < i < k. 



Proof. We give a constructive proof for this lemma. For a positive integer 
n = a ■ p e with gcd(a,p) = 1, we denote by u p (n) = £ the p-evaluation of n. 
Choose indices 1 < i < m and 1 < j < n so that 

v p (aij) = min{f p (o tK ) | 1 < i < m, 1 < k < n} 

By permuting the rows and columns of A, we may assume that ciij = an holds. 
Then an = a ■ p u p( ail > for some integer a with gcd(a,p) = 1. In particular, the 
integer a is invertible modulo pr. Multiplying the first row by the inverse a of 
a yields that the entries an, with 2 < i < m, are all divisible by aa\\ = p^ ail > 
(mod pr). Similarly, the entries aa\j, with 2 < j < n, are all divisible by aa\\. 
Therefore, we can find matrices L\ and R\ so that 



L x ARx 



I p Man) o • • • 0\ 


: A 

V o J 



holds. Continuing with the (m — 1) X (n — 1) sub-matrix A recursively would 
finally yield matrices L and R so that LAR = D is diagonal. A permutation of 
the rows and columns of D would then give the divisibility claimed above. □ 

The algorithm of Section U] now readily generalizes to an algorithm for solving 
the linear equation Ax = b (mod p e ): By Lemma [21 there exist matrices L and 
R so that 

LAR = diag(di, . . . ,dk,0, ■ ■ ■ ,0). 

Therefore the linear equation Ax = b (mod p ) translates to the systems Dy = 
Lb (mod p e ) and x = Ry (mod p ). It remains to check whether or not the 
diagonal system Dy = Lb (mod p ) has a solution and, if so, to determine 
all these solutions. This is a straightforward application of elementary number 
theory. 

Clearly for solving the linear system Ax = b (mod p e ) the explicit compu- 
tation of L is not necessary. Instead we can apply the row operations directly 
to the right-hand-side b which is cheaper in general. 

6.2 Solving linear equations using the block structure 

The algorithm of Section 16.11 generalizes to a method for solving a linear system 
of equations over arbitrary abelian p-groups as we will describe in the following. 

Let G and H be arbitrary abelian p-groups. Further let ip: G — > H be a 
homomorphism and let g € H be given. We decompose the groups G and H 
into direct products 

G = (Zp«i ) 91 x • • • x (Z p e e )« and H = (Zp-i ) hl x • • • x (Z p <* ) he 

with ei < ... < ei and <7j,/ij G N U {0}. Clearly, in the descending chains of 
subgroups 

G>p ei G> ... >p ee G = {0} and if > p ei G > . . . > p ee H = {0}. 



the factors p ei G/p ei+1 G and p ei H/p ei+1 H are homocyclic of rank gi + \ + • • • + gi 
and /ij+i + • • • + hi, respectively. Every homomorphism ip-.G^-H maps the 
subgroup p ei G to p ei H and therefore induces a homomorphism tp W : G/p e ' +1 G — > 
H/p ei+1 H Let ij, <5j, K{, and £j denote the natural homomorphisms so that the 
diagram 

G: 




G/p ei G 



^H/p e *H 



commutes. Since the quotients G/p ei G and H/p ei H are both homocyclic, the 
algorithm of Section 16.11 computes an x\ € G so that xf — b € ker k\ holds 
together with a generating set for some K\ C G with K^ 1 = ker ^W . In the 
remainder of this section, we lift these solutions recursively to a solution of 



;r T 



b. The following lemma lifts the special solution X{ 



Lemma 3 The solution x% lifts to a solution Xi + \ if and only if there exists 
X £ G such that X H € ker tpw with £1+1 = X{ — A and 



( x f _ M K »+i = x^Ki+i — y.+n' ( ' +1) ] 
Proof. Let A £ G be as above. Then 



(4) 



"i+l 



6)«*+i = (( Xi - A)* - 6) Kl+1 



UJ/tj + l 



^K; 



and hence xf +1 — b lifts the solution Xi- Assume that the solution X{ lifts to the 
solution Xj+i; that is, we both have xf, 1 — b E ker Kj + i and xf — b£ ker m. Recall 
that ker^j = p ei H > p ei+1 H = ker/tj+i. This yields that (xj+i — Xj)^ £ ker k, 
and hence 



o = (x i+1 - Xl y^ 

Therefore (xj+i — Xi) H € ker^W. 



(•?' 



i+l 



\t»¥> 



(0 



□ 



It remains to lift the kernel of the linear equation. Let Ki C G be a pre-image of 
ker y?W in G. Then a generating set for .fQ is easily obtained from the generating 
sets of keryjW and kertj. Suppose that Ki = (k±, . . . ,k() holds. Since we have 
that mpW = (fKi, it follows that 



It suffices to check whether or not xf 



b e (kf Ki+1 ,. . . , kp +l ) holds. Since 
j£<pk 1+1 ^ p e iH/p ei+1 H : the latter condition is equivalent to solve an equation 
over the homocyclic group p ei H/p ei+1 H. A solution to this latter system yields 
a lift of the solution Xj as described in Lemma [3j Furthermore, the following 
lemma outlines the lift of Ki to -JQ+i. 



Lemma 4 Let K- h < G be given so that K^ = ker y>w holds. Then it holds that 
K i+1 = ({k G K t | W K ^ = 0} U {p e *+i- e ' Jb | k G if<}). 

Proof. Let k £ Kibe so that fc^ re *+i = 0. Then, as </?Kj+i = ii + \Lp^ +1 \ it holds 
that fe ti + l¥ ' (l+1) = and therefore, k^ 1 G ker<p( i+1 ). Write Ap = p e *+i- e » and 
let fe G Ki be given. Then it follows that 

{Apk) H +^ +1) = (Apk)' pKi + 1 

and, since k^ G p 6i .ff = ker Ki, we also have that 

(Apjfc)* = ApW ep e ' +1 H = kerwi+i. 

This yields that L = ApKi is contained in -£Q+i and, as A;^ G kerKj + i < ker/tj, 
it follows that .JQ+i < i*Q. Hence, every element g G -fQ+i \ L can be written as 
g = ab with b G L and a G -fQ+i with aL ^ L. Then we get 

= n^+i^ (l+1) = (abY*^ 1 = df Ki+1 b VKl+1 

where W K%Jrl = 0. Hence the element a G Ki satisfies a ipKl+1 =0. □ 

Note that the generating set of K%+\ defined in Lemma3]may contain redundan- 
cies. These significantly affect the complexity of the algorithm in Section 16. fl 

6.3 On a number theoretical approach 

In this section we describe a number theoretical approach which avoids the 
redundancies introduced in the lifting of the kernels in Section 16.21 Our overall 
strategy for solving (TjQ) is an induction on the exponents fi, ■ ■ ■ , f m - More 
precisely, we solve the linear system Ax = b (mod p) over the finite field ¥ p 
with Gaussian elimination. Then we lift the obtained solutions recursively by 
applying the Hensel lemma. First we only consider an abelian p-group G and 
an endomorphism 92: G — )■ G. 

6.3.1 Solving endomorphic equations over finite p-groups 

Let G be a finite p-group and let <p: G — > G be an endomorphism. Further let 
b G G be given. We describe an algorithm for solving the linear equation x v = b 
or, equivalently, the linear system of modular equations 



a\\x\ + ■•■ + a lm x m = 61 (mod/ 1 ; 



(5) 



a m ixi + ••• + a mm x m = b m (mod/ m ). 

Note that the matrix A = (dy)i<i j< m in © satisfies the condition 

p e i~ e ™in{i,j} J a ,j for each 1 < i,j < m, (6) 

as it corresponds to the endomorphism (p; see also [8|. Denote the linear system 
of equations in ([5]) by Ax = b (mod [p ei , . . . ,p £m ]). 



Clearly, using Gaussian elimination, we can find all solutions to the linear 
system Ax = b (mod \p, . . . ,p\) efficiently Every solution x = (cci, . . . , x m ) to 
the system over the finite fields ¥ p has the form x = £o+£i^i + ' • -+£ r tr where, for 
each 1 < % < r, it holds that A£i = (mod [p, . . . ,p]) and ti G {0, ... ,p — 1}. 

The overall idea for solving ([5]) is to lift these solutions simultaneously by 
keeping the coefficients t±,...,t r in the finite field ¥ p . This yields that, if 
|ker</?| = p r , we obtain r independent homogeneous solutions £i,...,£ r . For 
£ G {1, . . . , e m } and k = min{i | £ < a}, we assume that all solutions to the 
linear system Ax = b (mod \p ei , . . . ,j? e&_1 ,j? , • . . ,p ]) are given by 

Z = 6) + tl£l + "- + tr&., (7) 

where ij 6 {0, ... ,p— 1} for each 1 < z < r. Applying the Hensel lemma, we lift 
the solutions in ([7]) to solutions to Ax = b (mod [p ei , . . . ,p 6fc - 1 ,p +1 , . . . ,p ]). 
For this purpose, we consider the modular equations 

afcixi + ••• + a fcm x m = 6 fc (mod/ +1 ) 

a-mixi + ■ ■ ■ + a mm x m = b m (mod / +1 ). 

By condition (j6]), the solutions »j (mod p ei ), for each 1 < i < k, do not need to 
be considered anymore. Similar to Section [67T], we can find a matrix L G Z mxm , 
which is invertible modulo p £+1 , so that A = LA is the matrix of an equivalent 
system of modular equations so that A = (aij)k<i< m ,i<j<m still satisfies 

p e i~ e tnin{i,j} I jj„ f or gach k < i < m and 1 < j < /c, 

and the sub-matrix {aij)k<ij< m is upper-triangular. Assume that all solutions 
to Ax = b (mod [p ei , . . . ,p ek ~ 1 ,p i , . . . ,p e ,p i+1 , . . . ,p e+1 ]) are given uniquely 
by 

Z = Co + Hi + ■■■+ t^r (8) 

with ti G {0, . . . ,p — 1}; that is, the solutions x m ,x m _i, . . . ,x n -\ modulo p 
are already lifted to solutions modulo p e+1 . We show how to lift the solutions 
x n = x n (ti, . . . , t r ) modulo p modulo p +1 . Consider the modular equation 

a/ixi H h aikXk + ai n x n -\ h ai m x m = bi (mod p e+l ) 

for some k < I < m. This equation can be considered a polynomial / G Z[x n ] 
with zeros modulo p^ given by x n = x n (ti, ■ ■ ■ ,t r ) for any choice of t\, . . . ,t r G 
{0, . . . ,p — 1}. More precisely, we define f(x n ) = f(x n ) — bi where 

f(x n ) = anxi H h a ik x k + ai n x n -\ \- ai m x m . 

The Hensel lemma now applies to the polynomial / and its zeros x n (ti, . . . , t r ): 
Let f'(x n ) = din denote the formal derivation of /. Using the linearity of the 
polynomial /, we obtain the following lemma which gives a condition for a 
unique lift of the zeros x n = x n (t±, . . . , t r ). 



Lemma 5 If f'{x n ) ^ (mod p) holds, then the solutions x n = x n (ti, . . . ,t r ) 
lift uniquely to x n + t ■ p e where t € {0, . . . ,p — 1} is given by 

i = -(a* n ) ^ 1 t— *H 1 7— *r (modp). (9) 

\ p l p z p l J 

Proof. The proof follows immediately from the Hensel lemma. □ 

Recall that a solution x n (ti, . . . , t r ) modulo p £ may not lift to a solution modulo 
p +1 . The following lemma gives a sufficient and necessary condition for such a 
lift. 

Lemma 6 If f'{x n ) = (mod p) holds, then the solutions x n = x n (ti, . . . ,t r ) 
lift if and only if the coefficients t\, . . . , t r € {0, . . . ,p — 1} satisfy 

- 1 7— h-\ 1 j-t r = (modp). (10) 

pi pi pi 

If this is the case, then the solutions x n = x n (t±, . . . ,t r ) lift to x n + t r +ip for 
each t r +i € {0, . . . ,p — 1}. 

Proof. The proof follows immediately from the Hensel lemma. □ 



The linearity of (1101) allows to eliminate a free variable t s , say, if this equa- 
tion is not trivially satisfied; that is, at least one coefficient does not vanish 
modulo p. More precisely, writing ji = f(£i)/p, for each 1 < i < r, and 
70 = (/(Co) - b)/p e , we may have that /(£,) ^ (mod p i+l ) but /(£ s +l) = 
• • • — f(£r) = (mod p +1 ). Then we may restrict 

t s = -77 1 (7o + 7i*i H h 7s-i*3-i) ( mod P)- (11) 

The following lemma determines an independent subset of solutions of ([HJ) which 
lift by Lemma [6j 

Lemma 7 ///(£,) ^ (mod p m ) fio/ds *Ie /(6+i) = • • • = /(&•) = 
(mod p), then the combinations of the elements 

{£,• — 7r 1 7»£s 1 fo?" each 1 < i < s — 1 
* s " " (12) 

£i, /or eac/i s + 1 < i < r. 

are p r ~ 1 solutions amongst f2j) which satisfy f{x n ) = (mod p +1 ). 

Proof. It is easy to see that for any element £j in (|12p it holds that f(x n ) = 
(mod p ). Thus it remains to prove that these elements are independent 
modulo [p ei , . . . ,p ek ,p e , . . . ,p i+1 ,p e+1 , . . . ,p +1 ]. For this purpose, assume that 
we are given ti, Ui € {0, . . . ,p — 1} such that 

fo + Yl **& = Co + 5Z Ui & ( mod [• • • ' p£+1 '■■•]) 



holds. Then, by construction, we have that 

Co + $^*i& - 7^(70 + 71*1 H h7s-i<s-i)6 

= Co + ^ u *& -7 s _1 (7o + 7i^i H 1- 7s-i«s-i)6 (mod [. ..,jp m ,...]) 

By Equation (|11|) . there exist <5, e G Z so that 

-7s" 1 (7o + 7i*i H 1- 7s-i*s-i) + Sp 



u s = -7 s _1 (7o + 7i^i H t-7 s _i'W s _i) + ep. 



(13) 



Thus we obtain that 

r r 

Co + Y ti & + 6 p^ s - £° + Y Ui & + e ^& ( mod [• • • >p* +1 » ■■■])• 



t=i t=i 

and, in particular, 

r r 

Co + ^ Uii + <Sp& = Co + Y u ^ + e ^ (mod [...,;/,...]). (14) 
i=i t=i 

If pCs vanishes modulo [. . . ,p e , . . .], the claim follows from the independence of 
Ci, • • • ,Cr m EJ Otherwise, as f(p£, s ) = (mod p e ) holds, the element p£ s can 
be written uniquely aspCs = «s+iCs+i + ' ' • + OrCr with a.% G {0, . . . ,p— 1}. Note 
that the the elements Ci, with s < % < r, have a weight strictly less than Cs! i- e - 
if if (C) = min{i | p*C = (mod [. . . ,p , . . .])} denotes the weight of C, then, by 
construction, we have that w{oti) < w(Cs)- Suppose that ^(oj) = w{^ s ) holds 
for some s < i < r. Then we obtain that 

as+i^- 1 Cs+i + • • • + a r p a ~ x Cr = p a Cs s (mod [. . . ,p e , . . .}) 

which contradicts the independence of £ s +i, . . . , Cr of weight p a if at least one on 
does not vanish. Thus, by induction on the weight, we obtain that Equation (11411 
implies that t\ = U\, . . . , t s = u s , and additionally 

ti + 5ai = Ui + eoti (mod p). 

By Equation (|13p . this yields that 5 = e and therefore, by the choice of the ij's 
and Uj's we finally obtain ij = Uj, for each 1 < i < s or s < i < r. □ 

We summarize the algorithm in 16.3.11 Note that we only used the fact that the 
linear system arises from an endomorphism of an abelian p- group, in form of the 
divisibility condition in (JU]). Clearly, the divisibility condition is always satisfied 
for any linear system of equations over a homocyclic p- group. Therefore the 
algorithm HENSELLlFTING also applies to the following problems: 

• Let G be an abelian p-group and let gi,...,g n G G. Decide whether or 
not b G (gi,... ,g n ) holds. 



10 



Algorithm HenselLifting ( A, b, p, e% < . . . < e m ) 
Determine the solutions to Ax = b (mod p). 
for £ G {1,... ,e m } do 

Echelonize the sub- matrix A modulo ]/. 

Let k(£) = min{i | £ < e^}. 

for n € {m, m — 1, . . . , &(■£)} do 

Let ctj n be the corner-entry corresponding to x n . 
if otj n ^ (mod p) then 

Evaluate Equation Q and determine t G {0, . . . ,p — 1} 
Lift the solution x n uniquely. 
else 

Consider the condition (|10p . 
if this cannot be satisfied then 

return fail; 
else 

Restrict to the solutions in Lemma [6j 
Lift these latter solutions. 
end if; 
end if; 
end for; 
end for; 

Algorithm 6.1: Solving endomorphic equations with the Hensel lemma 

• Let G be an abelian p-group and let <p: G m — > G n be a homomorphism. 
Then the algorithm HENSELLIFTING also solves the linear system of mod- 
ular equations Ax = b, as the divisibility condition is also satisfied in this 
case. 

• Let G be an abelian p-group with generators g\,...,g n . Further let 
ip:G — > G be an automorphism given by the images gf,...,g%. Then 
the algorithm HENSELLIFTING easily modifies so that it computes the 
inverse (p~ . 

The first problem yields that we can modify the algorithm in Section 16.21 for 
solving a linear system x^ = b for an arbitrary homomorphism ip: G ^ H of 
two abelian p-groups G and H. 

6.3.2 Solving linear equations over arbitrary finite p-groups 

Let G and H be arbitrary abelian p-groups and let ip: G — > H be a homomor- 
phism. The algorithm of Section 16.3.11 generalizes to an algorithm for solving 
the linear system x v = b for some b € H. This is a slight modification of the 
algorithm in Section \6. 21 
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